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D . Abstract: We explore the implications for neutrino masses and mixings within the 

minimal version of the supersymmetric large-extra-dimensions scenario (MSLED). 
' This model was proposed in hep-ph/0404135 to extract the phenomenological impli- 

. cations of the promising recent attempt (in hep-th/ 0304256) to address the cosmo- 

logical constant problem. Remarkably, we find that the simplest couplings between 
brane and bulk fermions within this approach can lead to a phenomenologically- viable 
pattern of neutrino masses and mixings that is also consistent with the supernova 
bounds which are usually the bane of extra-dimensional neutrino models. Under cer- 
tain circumstances the MSLED scenario can lead to a lepton mixing (PMNS) matrix 
close to the so-called bi-maximal or the tri-bimaximal forms (which are known to pro- 
vide a good description of the neutrino oscillation data). We discuss the implications 
of MSLED models for neutrino phenomenology. 
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1. Introduction 

The supersymmetric Large Extra Dimensional (SLED) scenario was proposed |l[, [^ in 
an effort to provide the long-sought understanding of why the vacuum energy should 
have incredibly small value — p^ac = -^^ with A ~ 10^^ eV — which is observed [0]. 
According to this proposal all observed particles reside on a 3-brane situated within 
an extra-dimensional world which has two large (1/r ~ 10^^ eV) dimensions^ (plus 

""^In our conventions the extra-dimensional volume is V2 = r^, and so (for a square 2-torus, for 
example) the Kaluza-Klein (KK) masses are of order m^^ ~ 27r/r. 
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possibly others which are much smaller, 1/r > 1 TeV), which can be possible 
provided the scale of gravity in the large extra dimensions is of order the TeV scale. 

The rationale for such a picture is that although non-gravitational (brane) par- 
ticles and interactions behave in the usual way, the gravitational response to their 
zero-point vacuum energy changes because gravity becomes 6-dimensional at energies 
above 10~^ eV. Although the ultimate success of the proposal is still under active 
study ^, 1^, its prospects for success remain promising. 

A particularly attractive feature of the SLED scenario is that it is very predic- 
tive, with many testable implications for cosmology and for tests of gravity on large 
and small scales. Better yet, it has observable implications for high-energy collider 
experiments, whose detailed implications have only recently begun to be explored 
1^. The predictiveness of the SLED picture is most pronounced in its minimal 
version, MSLED ||lOl, which assumes the particle content on our brane to be given 



by precisely the Standard Model particle content. 

Because of the similarity between the Kaluza-Klein (KK) scale in these models, 
27r/r ~ 10^^ eV, and the neutrino masses which are observed in oscillation experi- 
ments, it is natural to examine whether these models might also have observational 
implications for neutrino physics. We examine this issue here, and find a rich set of 
phenomena can arise within the neutrino sector. In particular, we argue that the 
observed pattern of neutrino masses and mixings can emerge through the mixing of 
the brane-bound neutrinos with the many massless 6D fermions which SLED mod- 
els predict must appear in the bulk, without running afoul of the bounds on sterile 
neutrinos. 

1.1 Neutrino Masses in MSLED 

The possibility that large extra dimensions might have implications for neutrino 
masses has been explored with many variations [^, [T^, |T^. The vast majority 



of these variations share the common feature that the neutrino mass eigenstates arise 
when the usual brane-bound neutrino flavours mix with various higher- dimensional 
neutral fermions which are assumed to reside within the bulk. It is generically 
assumed that these bulk fermions are massless in the higher-dimensional sense, and 
so from the 4-dimensional point of view only acquire masses (before mixing with the 
brane neutrinos) through the Kaluza-Klein mechanism. Various models are obtained 
by adjusting the strength of the brane-bulk couplings, the number of extra dimensions 
and the size of each dimension. 
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We follow a similar construction here, but find that the results from supersym- 
metric large extra dimensions are much more predictive, for the following reasons. 

• The massless field content in the bulk is dictated by 6-dimensional supersym- 
metry, and automatically contains numerous candidate neutral fermion fields 
which can mix with the brane neutrinos. Better yet, there is an understanding 
why the 6D mass of these bulk fermions must vanish. They vanish because 
they are related by supersymmetry to massless bosons, such as the graviton. 

• There is little freedom to alter the size and shape of the large extra dimensions, 
because the value of the observed Dark energy density dictates that there are 
two large dimensions, and the spacing of KK masses in both of these dimensions 
is not larger m^f^ ~ 2ii/r ~ 0.1 eV. 

• Like for other models with large extra dimensions, there is a natural exten- 
sion of lepton number to extra-dimensional fermions — chirality in the two 
large dimensions — which is also naturally broken at the scale mnn. This 
allows a technically natural understanding of why quantum corrections do not 
destabilize the observed neutrino mass pattern. 

• The back-reaction of the branes onto the bulk geometry is important in SLED 
models because it is part of the mechanism for suppressing the vacuum energy. 
In particular, this is how the bulk learns about supersymmetry breaking on the 
branes because the boundary conditions which the branes set up for the bulk 
fermions typically remove the massless fermionic Kaluza-Klein modes (like the 
gravitini) from the spectrum. This absence of bulk-fermion zero modes also 
removes the massless modes from the bulk fermions which mix with the brane 
neutrinos. This is important phenomenologically, since it is the degeneracy 
(in the absence of brane-bulk couplings) between these massless bulk states 
with the massless brane states which typically leads to large mixings between 
these states once brane-bulk couplings are turned on. In this way SLED mod- 
els sidestep one of the phenomenological difficulties of large-extra-dimensional 
neutrino models.^ 

Remarkably, we find that this restrictiveness can nonetheless very naturally lead to 
a phenomenologically successful picture of neutrino mixing. 

^Back-reaction issues are included in the more detailed microscopic construction of ref. UM- 
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1.2 Required Neutrino Properties 



Before describing the neutrino mass and mixing pattern to which SLED models lead, 
we first pause to summarize the main results of current neutrino phenomenology 
which must be explained. 



Active Neutrino Oscillations 



The observed pattern of neutrino oscillations are consistent with there being just the 3 
active neutrino species, whose participation in the charged-current weak interactions 
has the form 

igUa 



V2 



where 7^ is the projector onto left-handed spinors and i = 1, 2, 3 labels the 3 neutrino 
types, as a = 1,2,3 does for the leptons: {^1,^2,^3} = {e, /i,T}. Here the PMNS 
matrix [|16|, Uai, describes the amplitude with which the neutrino type 'i' reacts with 
the charged lepton 'a'. For three types of neutrinos it may be parameterized in terms 
of mixing angles and phases: U = V K, with 
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As is well known, assuming only three neutrinos to be relevant implies the split- 
ting between two of these neutrinos must be much smaller than their common split- 
ting from the third, and we follow common practice by choosing our labelling for 
the neutrino eigenstates such that the small splitting is between z/i and z/2. With 
this choice, the successful description of the oscillation lengths seen in solar- and 
atmospheric-neutrino experiments implies the three neutrino masses, mi^2,3; niust 
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satisfy 0, [18 1 

^ml^ = \ml-ml\ = (1.1 - 3.4) x 10^^ eV^ 

4 



Am| = \ml - ml\ = (5.4 - 9.4) x 10"^ eV^ . (1.4) 



This range for solar neutrinos corresponds to the Large-Mixing-Angle (LMA) solu- 
tion, which is the only one which accounts for the most recent results from the SNO 
and KamLAND oscillation measurements. In themselves neutrino oscillations do not 
fix the values of each of the masses separately, leaving open two possibilities. If the 
degenerate pair, z/i and z/2, are less massive than z/3, then the mass pattern is known 
as a 'normal' hierarchy, while if mi, 1712 > rris the hierarchy is 'inverted'. 

With this choice, the absence of the observation of short-distance oscillations of 
reactor-generated neutrinos implies the 3a limit 

|[/e3|' = 4< 0.067. (1.5) 

Because this implies C13 > 0.966, we have 1 > \Uei\^ + |t^e2p = c^s > 0.933 and so 
it is convenient to write \Uei\ = cuCrs ~ cos6'0 and \Ue2\ = suCis ~ sin6'0 when 
describing other oscillation measurements. Observations of solar-neutrino oscillations 
constrain S12, and imply 

0.67 < sin^ 2^0 < 0.93, (1.6) 

where, as above, 6q ~ 612. Oscillations as seen in atmospheric neutrinos are close to 
maximal, and satisfy 



0.85 < sin' 2^atm < 1 , (1.7) 



where 6'atm = 6*23- 



Suggestive Textures 

For later convenience it is worth remarking that a good first description for the 
observed values of the elements of the PMNS matrix is obtained by taking ^13 ^ 
and so Si3 ~ 0, Ci3 ~ 1. We may also approximate atmospheric- neutrino mixing 
as being maximal, in which case 623 ~ vr/4 and S23 ~ C23 ~ l/"\/2 = 0.707.., so 
sin' 2^atm = sin' 2^^23 = 1- With these choices the mixing matrix takes the form 

C12 S12 ^ 

-si2/v^ C12/V2 I/V2 . (1.8) 

^ S12/V2 -C12/V2 1/V2 J 



U 
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These choices are easy to obtain within specific models hke the extra- dimensional 
ones to be discussed shortly, because they follow as consequences of the assumption 



that the masses are invariant under a discrete Zi symmetry |]I9|, BO, |21[ which inter- 



changes the second and third generations. To see this, notice that the most general 
3x3 mass matrix which is invariant under this interchange may be written 



Z2 



mm m 
m' m + m m 
m' m m + m 



\ 



1.9) 



and this matrix always admits the eigenvector u oc (0, 1, —1)^ with eigenvalue m. As 
we shall see, imposing this symmetry on extra-dimensional models gives a matrix of 
this form, but with m = 0, leading to an inverted-hierarchy mass pattern for which 
z/3 is massless. 

In later sections we shall be led towards two particular choices of mixing matrices 
whose entries are very similar (up to signs) to the form of eq. ( p..8|) . For one case 

-l/VS : 



S12 



-0.577... and cu ^ ^/2/3 = 0.816... (and so sin 2612 



4^2 r2 



8/9 = 0.888...), which is called the tri-bimaximal-mixing form [22], and lies within 
the observed range allowed by oscillation measurements. In this case the PMNS 
matrix becomes 

/ -1/^3 ^ 

i/Vs 1/V2 . (1.10) 

A second example to which later sections lead is the bi-maximal mixing form, 
for which — S12 = C12 = l/\/2, and so for which the PMNS matrix becomes 



U 



( 1/^2 -1/^2 ^ 

1/2 1/2 l/y/2 
^-1/2 -1/2 1/v^y 



:i.ii) 



This form does not successfully describe the data, because the maximal value sin^ 26q = 
45^2^12 = 1 lies outside of the experimentally-allowed range. However it naturally 
arises if z/i and 1^2 are 'pseudo-Dirac', or carry an approximately-conserved lepton 
number |2^. In this case it is the perturbations which break the relevant lepton 
number which change S12 and bring the resulting PMNS matrix into agreement with 
observations. 
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Bounds on Individual Sterile Neutrinos 



Any 4D fermion which transforms as a singlet under the Standard Model gauge group 
can become a sterile neutrino if it mixes with any of the three active neutrinos. If 
such mixings exist then the PMNS matrix, Uai, for the charged- current interactions is 
no longer square, since it has 3 rows (a = 1, 2, 3, for each charged lepton) but 3 + 
columns {i = 1,...,3 + if the 3 active neutrinos are supplemented by A^ sterile 
counterparts). By virtue of these new mixing matrix elements sterile neutrinos can 
have observable signatures, none of which have been observed to date. This section 
summarizes the observational bounds which follow from these potential signatures. 

Many of the phenomenological constraints on sterile neutrinos are performed 
under the minimal assumption that only a single species of sterile neutrino exists 



{i.e. N = 1) [|^, |2^, |2^. The bounds which result typically rely on the absence 
of oscillations between this neutrino and the usual active ones, and so the strength 
of the allowed mixing can depend sensitively on what is assumed about the mass 
difference, Aml^, between the sterile and the relevant active neutrino state and 
about the active-sterile mixing parameter, e ~ |sin^^s| ^1.^ We summarize some 
of these bounds here, and return in later sections to describe the somewhat stronger 
bounds which arise for the higher- dimensional models of present interest. 



• Solar Neutrinos: The agreement between solar-model predictions and the ob- 
served charged-current and neutral-current solar-neutrino fluxes at SNO con- 
strain the amount of flux which can be lost to sterile states. For generic neutrino 
mass differences larger than 10~^^ eV^ these bounds constrain the sterile-active 
neutrino mixing angle to be smaller than e < 0.1, where the precise bound de- 
pends on the mass difference. This bound improves to e < 0.001 for mass 
splittings in a narrow range around Aml^ = 10~^ eV^, and to e < 0.01 for 
10~'^eV^ < Amf4 < 10~^ eV^, due to the presence in these cases of resonant 
active-sterile MSW oscillations within the Sun. 

• Reactor Neutrinos: The absence of neutrino flux disappearance from reactor- 
generated neutrinos constrains e < 0.1 for /S.m\^ > 10"^ eV^. 

^The presence of sterile neutrinos can also modify non- neutrino precision electroweak measure- 
ments, such as by modifying the muon lifetime from which the Fermi constant is inferred p7| , 
although these are not yet competitive with more direct observables. 
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Atmospheric Neutrinos: Active-sterile neutrino oscillations can be excluded for 
e < 0.2 for the range 10"^ < Amj^ < 10'^ eV^. 

Supernova Neutrinos: Excessively large active-sterile anti-neutrino conversions 
within supernovae can drain flux from the supernova neutrino signal, and too 
much of this can be excluded from the observed signal from SN1987a. For 
resonant z/g ~ '^s oscillations, this can constrain mixings down to e < 0.01 
for mass splittings, Am^^, which are larger than around 10 eV^. However 
because neutrino densities themselves play a significant role in these resonant 
oscillations, there is a feedback mechanism through which efficient resonant 
active-sterile oscillations can change the neutrino densities and thereby rapidly 



turn off the resonance 26, 28 



Nucleosynthesis: The dominant cosmological constraint for sterile neutrinos 
whose masses are smaller than 1 eV comes from the requirement that there not 
be too many light degrees of freedom contributing to the universal expansion 
during Big-Bang Nucleosynthesis (BBN). Even if the primordial abundance of 
sterile neutrinos is small, this constrains the strength of active-sterile oscilla- 
tions because these can cause too efficient production of sterile neutrinos from 
the equilibrated active neutrinos. This latter condition can constrain neutrino 
abundances down to mass differences of order 10~^ eV^, with the strongest 
constraints (coming for the largest mass differences) being of order e < 0.1. 

Cosmic Microwave Background: Measurements of CMB temperature fluctu- 
ations constrain the energy density which can be present in neutrinos at the 
epoch of recombination, leading to the constraint Quh"^ < 10~^. Since this con- 
strains the energy density present in neutrinos, the mixing angles to which this 
allows constraints are larger for heavier neutrinos, being sensitive to mixing 
angles as small as £ ~ 0.001 for sterile neutrino masses of order 10 eV. 



Bounds on Extra-Dimensional Sterile Neutrinos 

Since large-extra-dimensional models contain entire KK towers of sterile neutrinos, 

much more restrictive bounds are possible. For the purposes of discussing these 
bounds we anticipate the next section's results, and regard the sterile states to have 
masses ^ cg/r, and mixings ee with active neutrino flavours, where i labels the 
modes and q and eg are dimensionless parameters which depend on the details of the 
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compactification. As we shall see from explicit diagonalizations, typically ^ g/ci 
where g is a. constant which is independent of C.. 

The bounds in this case come in two forms. First, the presence of enormous 
numbers of sterile KK modes having masses ranging upwards from itikk ~ 27r/r, 
opens the possibility of there being active-sterile neutrino oscillations into specific 
sterile states for a great variety of oscillation lengths. This allows the use of the best 
of the bounds listed above for individual sterile neutrinos, provided only that the 
relevant mass differences are included amongst the allowed KK neutrino masses. 

The second kind of bound relies on the enormous available phase space which 
the extra dimensions make available, since this allows bounds to be obtained from 
the absence of processes which radiate energy into sterile neutrinos incoherently, 
rather than through coherent active-sterile neutrino oscillations. These incoherent 
neutrino bounds are in addition to the similar bounds which constrain the amount 
of incoherent emission which may be tolerated into other species of bulk particles, 
such as gravitons. A closely related class of constraints come from the cosmological 
limits on the total energy which can be tied up at present in sterile neutrinos. These 
bounds are particularly potent for extra-dimensional models, for which enormous 
numbers of massive and stable sterile neutrino states can be present. 

Bounds of these second type are normally used to completely exclude the pos- 
sibility of explaining neutrino masses in terms of more than a single large extra 
dimension, since each additional large dimension enormously increases the amount 
of phase space which is available in bulk neutrino states ||29|, ^ |3l|, ^ |3^ . The 
main novelty of our results in this paper is our ability to provide a phenomenologi- 
cally viable description of neutrino masses using two large extra dimensions without 
running afoul of these energy-loss bounds. 

Amongst the most restrictive new bounds obtained in this way are 

• Supernova Bounds: The large number of KK states makes energy-loss bounds 
from supernovae very strong. Each KK mode is radiated incoherently {i.e. by 
vacuum oscillations) with rate Fj ~ {m/Eyr^, where m ~ g/r is a measure 
of the mass term which mixes the active and sterile neutrinos (see the next 
section for precise definitions), E is the neutrino energy and Fj^ is a typical 
active-neutrino emission rate. Since the number of sterile neutrino states having 
mass smaller than E is of order (i?r/27r)" for n extra dimensions, an estimate 
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for the total incoherent emission rate into sterile states for n = 2 is 



r- = Vr 




V 



(1.12) 



Refs. ^ require this to be at most 10 * of the active-neutrino rate, F 



which leads to the conservative bound on the dimensionless bulk-brane mixing 
of order gjl-n < 10~^. 

We call this a conservative bound because it relies fairly strongly on the present- 
day understanding of supernova explosions, which are arguably relatively poorly 
understood given the inability of current computer models to successfully ex- 
plode a star.^ Another approach is only to require that the rate of sterile- 
neutrino radiation not be more than of the same order as the active-neutrino 
signal, in which case only a comparatively weak constraint on g is obtained. 

For supernovae, another potentially very dangerous mechanism for energy loss 
into sterile neutrinos is possible because the large number of neutrino states 
allows a succession of potential resonant oscillations which can depress the 
active- neutrino survival rate as it passes through the supernova environment. 
For small mixing each resonance is narrow and the survival probability can 
be understood as the product of that for each KK mode. Even though each 
KK mode mixes with an amplitude which is suppressed by the mode's mass, 
the resonance for each has essentially the same adiabaticity parameter, leading 
to a prohibitively large conversion rate even for extremely small brane-bulk 
mixing parameter g. Happily, as pointed out in refs. this constraint is 

too naive since the feedback of such oscillations onto the supernova neutrino 
densities is likely to quickly reduce the rate of energy loss in this channel to 
the vacuum-oscillation rate described above. 

• Late-Epoch Cosmology: Large-extra-dimensional models face well-known dan- 
gers with cosmology since the various KK states for the many bulk fields can 
cause problems with Big Bang Nucleosynthesis or with over-closure of the uni- 
verse if they are too abundant |3^. These bounds can be avoided if the KK 
modes do not appreciably decay into photons and if they and their decay prod- 
ucts are not too abundant at the BBN epoch. Although it remains a challenge 
to obtain a pre-BBN cosmology which gives sufficiently few relic KK modes 
^We thank Marco Cirelli and George Fuller for helpful conversations on this point. 
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in models with large extra dimensions, we put this issue aside here since its 
proper understanding must also await a hitherto missing study of the time- 
dependence of 6D cosmological background solutions. We here therefore follow 
the usual practice and assume there to be an acceptable number of KK modes 
at BBN, and focus specifically on the neutrino-related constraints which arise 
due to the possibility of populating KK sterile neutrinos through active-sterile 
oscillations from the known equilibrium abundance of active neutrinos. 

Constraints on the present-day energy density in neutrinos preclude there be- 
ing a significant number of stable massive sterile neutrino states having masses 
greater than of order 1 eV and having an appreciable mixing with active neu- 
trinos. Mixings are dangerous in this context because they can keep sterile 
neutrinos in equilibrium down to low temperatures, leading to too large late- 
time abundances. An efficient way to evade these bounds is to require all of the 
sterile neutrinos to thermally decouple at a temperature. To, above the QCD 
phase transition so their abundance can be diluted by the entropy release which 
occurs during this transition. This requires To > 1 GeV, where for each sterile 
neutrino species Tu is related to the active-sterile mixing parameter e by the 
condition 

^^G|r^-^, (1-13) 

where Gf ~ 10~^ GeV~^ is the Fermi coupling and Mp ~ 10^^ GeV is the 
Planck mass. The condition Td ~ 1/ {e'^G'j.Mp)^/^ > 1 GeV leads to the con- 
straint e <! 10~^. In later sections we find £i ~ g/ci, where g measures the 
dimensionless brane-bulk coupling and Q/r > 27r/r is the mass of the KK mode 
in question. We see that the constraint is strongest for the lightest modes, and 
for modes with q ~ 27r this constraint implies g < 2tt x 10"''. 

We see that both of the above constraints may be satisfied provided that the 
brane-bulk couplings satisfy g <, 27r x 10~^, and for this reason we have the regime 
of small g in mind when we describe the neutrino mixings in more detail in the next 
section. 

2. Neutrinos in MSLED 

In MSLED neutrino masses can arise through mixing between the usual three brane- 
bound neutrinos, Ua, and fermions in the bulk. In order to describe this mixing we 
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first digress to summarize our conventions for 6D spinors. 



2.1 6D Spinors 

In 6 dimensions tlie Dirac matrices are 8x8 matrices, wliicli we take to fiave tlie 
following form: Fm — {r^, ^m}, where /i — 0,1,2,3 and m = 4, 5 and 

r;x = 7/i®^2 and r„ = 75 8) • (2.1) 
Here denotes the usual 4x4 Dirac matrices, while 



r4 = . and rg = . (2.2) 





In addition we define the usual 4- and 2-dimensional chirality projectors = ^(1 + 
75), 7ii = ^(1 - 75) and r± = |(1 ± rg), with 




^3=1 1 ■ (2.3) 



In terms of these matrices, 6D chiral spinors are distinguished by the eigenvalues of 
the matrix Tj = ® T3. 

With these choices 6D spinors carry 4D and 2D indices, {u, a), where u = 1, 
and a = 1, 2. It is useful to take the basis for these spinors to have definite chirality, 
since a 6D spinor with positive chirality, F^ip — has 4 complex components 
and decomposes into two 4D states for which the eigenvalues of the pair (75, T3) are 
(A, s) = (+1, +1) and (A, s) = (-1, -1). 

Once compactified, we see that each 6D fermion reduces to a Kaluza-Klein tower 
of 4D fermions. For example a collection of 6D = +1 Weyl fermions has a 4D 
expansion of the form 

N'{x, y)^lj2J2 nlix) , (2.4) 

where i labels the KK modes, with ues{y) being 2-component spinors which are 
chosen to satisfy T3 u^g — s u^g- Consequently (because T^N^ — +N^) the 4D spinors 
nig have 4D chirality 7571^5 = sn^g (notice that in the above all spinor indices are 
suppressed). An overall factor of the size of the extra dimensions, 1/r, has been 
extracted from the extra-dimensional fermion wave- functions, Uis{y), so that these 
can satisfy r-independent orthonormality relations. 
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2.2 Bulk-Brane Couplings 

The lowest- dimension bulk-brane interaction between left-handed lepton fields on the 
branes and a collection of = +1 6D fermions, A^^, has the form 



int,N 



J d\ [\aULlNi^)H + C.C.] . (2.5) 



Here the Xaia are coupling constants having dimensions of inverse mass, which we 
expect to be roughly of order in size, where Mg ~ 10 TeV [|1^ is the gravity 
scale in the bulk. In this expression the SUl{2) gauge index is suppressed, while the 
index a = 1,2,3 labels fermion generations. Similarly u = 1,...,4 and a = 1,2 are 
the 4D and 2D spinor indices, which we here temporarily re-instate. 

It is convenient to write the couplings in terms of definite 2D chirality, with 

with T3 M^^-* = su''^\ If both A^^^ and x[^j ^ are nonzero for some choice for / and a then 
these brane-bulk couplings violate the 2D local Lorentz transformations generated 

by 

■J = l r , r^] = '-[10 [n, n]) = i (/ ® rs) , (2.7) 

which act on the 2D spinor index a. On the other hand, if all of the A^^"* (or X^^j^) 
should vanish, then the coupling preserves this 0(2) sjTumetry provided that the 
brane fermions are also taken to be charged under it. Any such 0(2) symmetry 
would be naturally interpreted as a lepton-number symmetry for the neutrino mass 
matrix, since this is the only possible way such a symmetry can be extended to an 
invariance of the rest of the Standard Model brane couplings given the assumption 
of minimal particle content. 



Models with bulk-brane neutrino mixing have been studied in detail in ref. [11 



T^, O, M , but this earlier round of model building differs from the present picture 



in several important ways. 

• Brane Back- Reaction: First, these earlier analyses ignore the brane back- 
reaction onto the bulk. However, we know quite generally in MSLED that the back- 
reaction of the brane (or branes) onto the bulk geometry induces boundary conditions 
for the bulk fermions at the brane position, and these remove the otherwise-massless 
Kaluza-Klein level [|l|. (Indeed this is how the bulk modes in MSLED 'see' that 
supersymmetry is broken by the branes.) 
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• Approximate Lepton Number Conservation: We have seen that the brane- 
bulk mixing relates lepton number on the brane to local 0(2) Lorentz transformations 
in the bulk. Because of this connection it is guaranteed to be a symmetry of the inter- 
actions of the 6D supergravity in the bulk. In general this symmetry is spontaneously 
broken by the background geometry, since extra-dimensional Lorentz rotations are 
always broken by the zweibein of the internal two dimensions, e^". For some ge- 
ometries (such as a 2-sphere) it can happen that the transformation of under 
such Lorentz transformations can be compensated by performing a diffeomorphism 
in the extra dimensions, provided that the extra-dimensional geometry admits a ro- 
tational isometry. In such cases the background breaks the product of local Lorentz 
transformation and diffeomorphism symmetries down to the diagonal rigid rotation 
which preserves e^", in an extra-dimensional analogue of the spin-orbit coupling. 
For geometries without rotational isometries (such as 2-tori) no such compensation 
is possible and the symmetry is broken completely. This leads to a naturally very 
small bulk lepton- number breaking scale, of order the KK mass scale, 771^^ ~ 0.1 
eV. Because of this low symmetry-breaking scale, it is technically natural to have 
the scale of lepton-number violation in the bulk-brane couplings also be as small as 

• Bulk SUSY: In previous studies the properties of the bulk fermions were usu- 
ally adopted for the convenience of neutrino phenomenology, for lack of a theory of 
the origin of the bulk fermions, . In MSLED, however, supersymmetry strongly 
constrains the properties of the bulk. There are, for instance, fermions which are 
universal to 6D supergravity coming from the supergravity multiplet itself, which 
typically contain the dilatini, and/or the extra-dimensional components of the 
gravitino, ■0^^,, m — 4,5. If the fermions, , in the bulk-brane mixing terms are 
linear combinations of these fields, their bulk interactions are strongly constrained by 
the fact that they are related by supersymmetry to the extra-dimensional graviton. 
In particular, such fermions cannot have any explicit 6D mass terms at all. (Super- 
symmetry similarly constrains the properties of other bulk fermions, which appear 
in 6D matter multiplcts (hyperini or gaugini), since such fermions are also typically 
tied by supersymmetry to massless bosons in the bulk. 

2.3 The 4D Neutrino Mass Matrix 

We now write down the neutrino bulk-brane mixing which results in a basis of fields 
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for which we assume the charged-lepton masses to aheady be diagonal. Because 
no bulk 6D mass terms are allowed by supersymmetry, the complete 4D neutrino 
mass matrix consists of eq. ( |2.5| ) plus the Kaluza-Klein masses coming from the 
extra-dimensional kinetic terms of the bulk fermions. The result is 



£1 



{yalLn[+) + -^^{yalLn\_) + — {n[^-iLn\_) + c.c 



(2.8) 

where v = 246 GeV is the Standard Model Higgs v. e.v., and c| 7^ is an 0{2tt) dimen- 
sionless number whose value depends on the detailed shape of the extra dimensions 
and which controls the form of the Kaluza Klein masses through mj = c\/r. For 
instance for compactifications on a square torus (without brane back-reaction) £ is 
a pair of integers, £ = (fci, ^2), and q = Ck^k2 = 27r(A;i + ik2). Here also h denotes 
the 4D left-handed spinor which is the conjugate to the right-handed spinor n. With 
(2,0) — or higher — 6D supersymmetry in mind, we do not write an s = ± de- 
pendence for c| because we assume both 2D chiralities to share the same 6D kinetic 
terms and boundary conditions in the bulk. 

Allowing the indices to run over the ranges a = 1, 2, 3, s = ± and / = 1, . . . ,n, 
the resulting (left-handed) neutrino mass matrix becomes 









sSr' sir' si? sS' si;' 

921 921 922 922 923 923 

9S' 9ii' 9iP 9^' 9^ 9i^ 








c- 



(+) (+) (+) 

^11 921 931 

9ii' 9ii^ 9ti^ 4=' 



9l2^ 9^,2 932 
9l2 922 932 








1=1 

I 











J=2 








J=2 










V 










(2.9) 



where g[^] = X[^jV. For each i the matrix of diagonal 2x2 blocks with elements c| 
is a 2n X 2n matrix corresponding to the n bulk fermions, assuming these all share 
the same boundary conditions (and so also the same KK levels i). We do not require 
detailed expresions for the c|'s, but because of the brane back- reaction we do take 
all of the 4 7^ 0. 

In general this matrix is complex and symmetric, and so may always be diago- 
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nalized by finding an appropriate unitary matrix, 11,^, such that 



Hi 

/i2 

/is 



V 



(2.10) 



where the superscript 'T' denotes transposition [not hermitian conjugation). The 
required Ui, is 

= V,K (2.11) 

where is the unitary matrix which diagonahzes the hermitian matrix Mj^Mjy — 
but in the usual fashion of a similarity transformation: Vj (MlMt,) = diag(/i^, . . . ) 
— and 7^ is a diagonal matrix of phases which can be chosen to ensure that the /ij 
are non-negative. If the matrix M^, should be real — as we assume for simplicity 
below — then can be chosen more simply to be the usual orthogonal matrix which 
diagonalizes it. (This diagonalization is performed explicitly for various choices for 
the mass matrix below.) 

Given such a diagonalization, the three species of 3-brane neutrino states may 
be written as a linear combination of the mass eigenstates according to 



(2.12) 



where the non-square matrix U consists of the first 3 rows of the matrix 1/^- Putting 
this into the charged-current interaction, eq. ( |1 . Ij ) , shows that it is this matrix, 
U, which plays the role of the PMNS matrix in extra-dimensional models, and in 
particular the absence of evidence for sterile neutrinos requires its left-most 3x3 block 
to agree with the 3x3 PMNS matrix described above, which is inferred from neutrino- 
oscillation data. Similarly, its 4th and higher columns must all be sufficiently small 
not to conflict with the various bounds on the existence of sterile neutrinos. 















= u 








z/4 




['■■) 



3. Explicit Models 

In order to proceed further it is necessary to make some choices for the brane-bulk 
couplings X[^j, and for the shape of the internal dimensions (which determines the 
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coefficients c|). In this section we consider the simplest choices in some detail, guided 
by considerations of technical naturalness. For simplicity we assume only a single 
bulk field couples to the brane, and so henceforth suppress the bulk index J, and so 

3.1 A Toy Test Case 

Before laying out the couplings which we expect MSLED to provide and which furnish 
acceptable phenomenology, it is instructive to pause here to examine a simple choice 
for couplings which does not work, in order to identify better the issues which must 
be addressed by a successful choice. For this purpose consider the choice where all of 
the couplings are equal, regardless of flavour and chirality: Aa = A for all a = 1, 2, 3 



and s = ±. As is easily verifled, with this choice the mass matrix, eq. ( |2.9| ), leads to 
the following neutrino masses and mass eigenstates: 

• Massless States: Two massless states are given by 

z/o = ^(-2,l,l,0,---)'^ and z/^ = (o, 1, -1, 0, ■ ■ ■ . (3.1) 

Because the nonzero entries only appear in the flrst three positions, these massless 
states do not mix with the inflnite tower of bulk modes. 

• Massive States — Part I: For each £ there is a neutrino eigenstate having mass 
/i = \m(\ = \ci\/r, given by 

z/fc = -^(0,0,0,---,0,l,-l,0,--y, (3.2) 

where k = 1,2,3,... and the two nonzero values appear in the {2k + 2)'th and 
{2k + 3)'th entries of the column vector. These states clearly have the same KK 
masses as would have been present without the brane-bulk mixing term, and since 
the flrst three entries of the eigenstate vanish these states do not mix with the brane 
neutrinos. 

• Massive States — Part II: All of the other eigenstates also have masses which are 
nonzero and depend nontrivially on the coupling g = Xv. These states have masses, 
fii = \zi\/r, where zi are deflned as the roots of the following equation 

- = (3.3) 

The corresponding eigenstates are given by: 

/ 1 1 1 1 1 1 1 \^ , . 

V ^g ^g ci=x - zi q=i - zi q=2 - zi q=2 - zi J 
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with normalization constant A*";. 

It is instructive to ask what happens to these states in the hmit g ^ 0. As 
may be seen, for instance by solving this equation graphically, all but one of the /i;'s 
approaches the KK masses, \mi\ = \ce\/r in this limit. The one state whose mass 
does not approach \me\ as g ^ has a mass, /io, which is bounded above by 

/io < — \<S\ , (3.5) 
r 

where 



e ^ 

Evaluating /iq perturbatively in g shows that actually equals this upper bound up 
to corrections which are 0(^g^\ 

Although this mass nominally vanishes q& g 0, what is interesting is that 
the sum appearing in /io generically diverges in the UV, with a divergence linear 
in the cutoff L ~ Mgjm ~ IS/lgT I (Stt) ~ 10"'^^ for two extra dimensions — where 
the numerical estimate takes Mg ~ 10 TeV and TTijfjf ~ 0.1 eV. Keeping track of 
factors of 27r, and recalling q ~ 2n, we take S ^ ^ {2'kL) = L. This illustrates 
how extra-dimensional models can complicate the problem of identifying masses and 
mixings in the small-g limit, because the divergences associated with sums over the 
very large number of KK modes (like S) can overpower small pre-multiplying factors 
(like g'^), even for the very small values g'^/{2'K) < 10^^ of interest in what follows. 

The divergence appearing in /xq may be absorbed as a renormalization of a 
counter-term localized on the brane, which has the form of the usual Standard Model 
dimension-5 neutrino mass: 

^brane = hat J d^X {LaH){UH) . (3.7) 

The required renormalization is of order 5h ~ X^Mg ~ ^/Mg. Since the neutrino 
masses induced by such an interaction are of order hv"^ and the natural scale for 
the hab is of order 1/Mg, the existence of such a term underlines the fact that any 
successful explanation of neutrino masses using brane-bulk couplings within MSLED 
must also explain the absence of this pure brane term, which can produce much too 
large a mass [0 . 



A second type of divergent sum which we shall also encounter is g V, where 



Because this diverges logarithmically for two extra dimensions — P ~ -^^^{271 In L) ~ 
2^ In^Mg/mi^K) ~ 32/(27r) — the combination g'^V can be large provided g'^/{2n) > 
0.03. In practice, our interest in the regime g < 10""^ allows us to neglect such 
sums, but we nevertheless provide in an appendix expressions for the eigenvalues 
and eigenvectors which do not assume the quantity g^V to be small. 

3.2 Perturbative Eigenvalues and Eigenvectors 

As the above example shows, agreement with the neutrino data requires more com- 
plicated bulk-brane couplings than those which were entertained for the Toy Model. 
Before turning to specific proposals for more phenomenologically successful brane- 
bulk couplings, we pause here to record general expressions for the eigenvalues and 
eigenvectors for more complicated neutrino mass matrices. Keeping in mind the 
strong bounds on active-sterile mixing, and in order to keep the analysis simple, 
we solve the eigenvalue and eigenvector problem to leading nontrivial order in the 
brane-bulk couplings, after first stating some exact definitions. (We also provide a 
more detailed calculations for some phenomenologically interesting choices for the 
mass matrix in the appendix.) 

Our interest is in diagonalizing the mass matrix of eq. ( |2.9|) , in the special case 
of a single bulk field but with couplings X'^^ which are otherwise arbitrary. To this 
end we write the mass eigenstate corresponding to the lowest three mass eigenvalues, 
fJ-i = fJ'0,± = Zi/r = 0{g^/r) as follows: 

Ui = Ni(nieii,nie2i,nie3i,xii,yii, . . . ,xei,yei, . . , (3.9) 

where Ni is an overall normalization constant and rii is chosen to ensure that the Cai 
satisfy e^j + e^j + = 1. The normalization constant is given by A^j = [n^ + Ef]"-*^/^ 
where Sj denotes the sum 

^^ = J2{xl + yl). (3.10) 

e 

In analyzing the physical implications of any such an eigenstate it is convenient 
to define the active-sterile mixing angle by Csi = cos6si = Ni rii, so that the first three 
entries for Ui become c^j Cj, with Cj ■ Cj = 1. In terms of Sj this is may be written in 
the useful equivalent form 

J]. 

tan^,i = — . (3.11) 
rii 

Notice that the unitarity of the full matrix Uai, together with the unitarity of the 3x3 
matrix whose elements are Cai, implies that this angle is the one which is relevant to 
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the total incoherent energy loss into sterile neutrinos from a charged- current reaction 
involving lepton flavour ia'- 

Ra = J2 = 1 - E = 1 - = . (3.12) 

j>3 1=1 i=l i=l 

Writing, as before, = X^a^^v, it is fairly simple to compute the mass eigenval- 
ues and eigenstates to second order in the small quantities Within perturbation 
theory, the leading-order masses, /ij, and the components, Cai, for the lightest three 
states are determined by diagonalizing the following perturbative correction to the 
mass-matrix within the degenerate 3-dimensional massless subspace: 

f^a, = -^{9^-'^9l-^ + 9l^^9i-^)- (3.13) 

This matrix has rank 2, which ensures that one of the three light neutrinos remains 
massless to this order. Once the components Cai, are obtained in this way, the leading 
expressions for the components Xa and are similarly given by the perturbative 
expressions 

and ya = , (3.14) 



where Cj ■ = X]a=i ^ai 9c^^ ■ Similarly, the quantity rij = 1 up to corrections which 
are 2nd order in the 5'*^^-'. To leading order the active-sterile mixing angle, 9si, for 
the three neutrino species whose mass vanishes as g ^ then becomes 

tan^ ^ = 5^(4 + yl^i = [(el ■ g^^^f + (e. ■ g^'^f] V , (3.15) 
where V is as in eq. (|3.8|) . 



3.3 Approximate Symmetries I 

With these insights in mind, we next explore a more realistic model of neutrino 
masses within the SLED framework. We start in this section with a model having 
a minimal flavour content, which can satisfy all of the required phenomeno logical 
constraints but one. Even though it does not succeed in the end, we describe this 
model in detail here for several reasons. First, we do so because it comes so close and 
illustrates well the tension between the various constraints which makes constructing 
a successful model difficult. Second, it is useful to describe in detail the symmetry 
issues which the model raises, since these also play an important role in the successful 
model we describe in the next section. 
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The simplest and most natural choice to make for the brane-bulk couplings at 
the TeV scale is to suppose them to be both flavour-independent and Lepton-number 
conserving. Lepton number conservation is necessary because it precludes the ap- 
pearance of a brane neutrino mass like eq. ( |3.7| ). Flavour-independence is a natural 
expectation since couplings to the bulk are gravitational, and so the microscopic 
physics which underlies them is unlikely to care about the details of the low-energy 
flavour physics on the brane. 

These symmetries suggest we take the couplings at the TeV scale to be 

AW = and Ai-) = A, (3.16) 

for all brane flavours, a, which preserves the Standard-Model lepton number, L = 
L^ + L^ + Lt, in addition to the extra-dimensional chirality, J. With this choice there 
are precisely 3 massless neutrino eigenstates as well as a tower of KK modes whose 
Dirac masses are nonzero. The masslessness of the three lightest states is protected 
by the lepton number invariance. Notice that because these states are degenerate in 
the leading approximation, their mutual mixings need not be small once this mass 
matrix is perturbed to obtain more accurately the neutrino masses and eigenstates. 
We do not have more than 3 such degenerate states — unlike previous analyses — 
because of our inclusion of the brane back-reaction, which removes any KK zero 
modes. 

Now we imagine renormalizing these couplings down to the lower energies rel- 
evant to neutrino physics, and we expect that neither the flavour nor the lepton- 
number invariances are likely to survive this process as an exact symmetry of the 
low-energy theory. Flavour symmetry is broken at the very least because the Stan- 
dard Model brane couplings themselves distinguish amongst flavours. Similarly, we 
have seen that within the bulk the role of lepton number is carried by local Lorentz 
rotations, J, in the extra dimensions, and these are naturally broken by the back- 
ground extra-dimensional geometry, such as by the background zweibein, 6^"-^ The 
symmetry-breaking scale for this symmetry is then naturally of order the KK scale, 

^Notice that Cm" can remain invariant if the effects of the local Lorentz transformation acting 
on a is cancelled by performing a compensating diffeomorphism on the index m. The diagonal 
transformation which preserves the zweibein can then be an unbroken symmetry of the full theory 
if the compensating diffeomorphism is also an isometry of the metric — em"" Sn'' Sab- For 
instance, it can be unbroken if the internal metric is rotationally-invariant (like for the 2-sphere) 
but would be broken if not (like for the 2-torus). 
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ttlkk ~ 271 /r. Once such corrections are included into the brane-bulk couphngs we 
expect eq. ( p.l6| ) to be replaced by 

^i+) = ei+) and ^(-) = ^ + 4') , (3.17) 

where as before g = Xv, and the ei"^^ are of order mK^/Mj ~ Kitli^j^ / Mg ~ 
27C K/{Mgr) ^ (7 ~ "^/^gi where Mj denotes the mass scale which communicates 
the symmetry-breaking scale to the brane sector, and K = Mg/Mj is a dimensionless 
constant which we choose to obtain an acceptable neutrino phenomenology.^ Our 
task now is to infer the neutrino masses and mixing which are implied by such a 
coupling choice. 

To second order in g^"^ the effective mass matrix, eq. ( p.l3|) , which determines 
the mass shift of the three massless unperturbed states becomes 

9^ (A+) , .(+) 



(3.18) 

where the approximate equality drops terms which are of order in comparison with 
those which are order eg. 

This matrix has rank 2, which ensures that one of the three light neutrinos re- 
mains massless to this order. The other two neutrinos acquire masses which are gener- 
ically of order eg/r times the divergent sum S. Keeping in mind that e ~ Kitikk / Mg 
and the sum diverges linearly in the cutoff, which is of order L ~ Mg/mj^i^, we see 
that the nonzero masses which are generated are of order egL/r ~ gK/r. Because 
of the protection of the approximate lepton number symmetry — i.e. the small size 
of e^'' — these are the right order of magnitude to describe the observed neutrino 
masses.^ 



Obtaining the precise form for the masses and mixings of the lightest three 
states requires making more specific choices for the ei^^'s, and guided by neutrino- 
oscillation phenomenology described in earlier sections we choose these to preserve 



^As discussed in ref. |0, an explicit microscopic realization of the SLED picture is likely to 
involve a number of mass scales smaller than the 6D Planck constant, Mg, such as the string scale, 
Ms, and the KK scale for any dimensions in addition to the 6 discussed here. 

''To the extent that bounds on sterile neutrinos require us to take g < 10~^ we must imagine 
K ~ 10^. An explanation for why K should be so large is a requirement for any explicit SLED 
compactification which hopes to describe neutrino phenomenology. We here content ourselves with 
a phenomenological approach, showing that neutrino phenomenology need not conflict with sterile 
neutrino bounds in SLED models. We leave for future work the explicit calculation of symmetry- 
breaking effects like K in an explicit model. 
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the Z2 symmetry which interchanges the 2nd and 3rd generations. In particular we 
take 



.(+) 



ei,e2,e2 



in which case the one massless eigenvector is the standard one, 

1 / xT 



0,1, -1,0,--- 



(3.19) 



(3.20) 



for which X£o = ym = sin^^o = 0. The two nonzero eigenvalues which vanish as 
g become 



(ei + 262) ± ((ei + 262)2 + 2(6i - 62) 



1/2- 



(3.21) 



where V is as defined in eq. ( |3.8| ). Notice the n± are 0{eigS/r), as required. The 
corresponding neutrino eigenstates have 



L=(a±,l,l) 



where 



a/2 + u;_t 

(61 - 262) ± ((61 + 262)2 + 2(61 - 62) 



(3.22) 



1/2 



It follows that is given by 



tan^ Q„ 



ei + £2 



(a± + 2)- 
ai + 2 



g'v. 



(3.23) 



(3.24) 



where the approximate equality uses ei <^ g. 

Notice that if 61 = 62 = 6, then a+ = 1 and a_ = —2, and so tan^ 9s+ = ^g^V 
and tan^ 9s- = 0. This leads to the following PMNS matrix for the lightest three 
neutrino states 

^-^/2/3 Cs+/V3 ^ 

1/VQ Cs+/VS I/V2 , (3.25) 
^ 1/VQ Cs+/VS 1/V2j 

whose matrix elements agree (in the limit 6s+ 0) in magnitude with the tri- 
bimaximal mixing matrix described earlier. For these same choices the neutrino 



mass values of eq. (|3.21| ) become /i_ = and /i+ = 6egS/r. 

The problem with this model is that it does not appear to have a parameter 
range for which the neutrino mass spectrum, mixing angles and sterile-active mixing 
angle are all acceptable for the same choice of parameters. To see this notice that 
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if we imagine taking g sufficiently small to minimize the active-sterile mixing, then 
we expect the observed neutrino oscillations to be obtained from the mixings of the 
three states whose masses vanish in the limit g 0. Since our conventions assume 
that for three-neutrino mixing it is z/i and 1^2 which dominantly participate in solar- 
neutrino oscillations, we see that the observed hierarchy lAmU <^ |Am^^j^| then 
requires us to take — /i-)^ ^ |yU+yU_|, and so |ei + 2e2| ^ 2|ei — 62!. Notice 
that this condition is precisely the opposite to the choice ei = €2 which led to the 
successful PMNS mixing matrix, above. 

Although we have not found a choice for the which produces acceptable masses 
and mixings using the perturbative analysis just described, we can do so if we make 
g large enough that the quantity g^V is not small. In this case the more exact 
diagonalization of the appendix shows that eq. ( p.21|) generalizes to fj,± = \z±\/r, 
with 



z±= gS 



(ei + 262) ± ((ei + 2e2)2 + 2(ei - 62)^(1 + 3g^V) 

1 + 3g^V 



1/2' 



(3.26) 



and so the condition for obtaining a small ui — 1/2 mass splitting becomes 



(ei + 262)' <2(ei- 62)^(1 + 3/P). 



(3.27) 



When g'^V > 1 the quantities which control the mixing angles become 



2[ig^ + eie2)ziV{z^) + ^(ei + e2)5(zi)] 2[g^z,V{z,) + g{ei + t2)S{zi)] 



Zi - {9'^ + el)ziV{zi) - 2geiS{zi 



Zi - g^z.{P{zi) - 2geiS{zi 



(3.28) 



and 



yei 



g{ai + 2)ci + {aiEi + 2e2)zi g{ai + 2)ci 



- 4 



Za Co 



- cj 



g{ai + 2)zi + (a^ei + 2e2)Q g{ai + 2)zi 



2 2 ' 

Z: Cn 



(3.29) 



where the approximate equality for yu applies for all £ if Zi is not too small, but is 
not valid for i too close to the cutoff if Zi ~ 0(1). Here, as before, \zi\ = ^ir, and 
the functions S{z) and Vi^z) are defined by the sums 



S(z) 



z"^ - 



and V{z) = J2 -TZ 



(3.30) 
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(from which we see 5(0) = —S and 'P(O) = —V, and so in particular eigS{z) ~ 0(1)). 
Finally, assuming the sum over x^^ + yl^ to be dominated by its UV divergent part, 
we have 

g^V. (3.31) 

Prom these expressions we see that /x+ and can be made sufficiently degenerate 
without giving up an acceptable PMNS mixing matrix, e^j, provided wc take g'^V to 
be sufficiently large. (Recall that P ~ 32 and so g^V can be taken to be large while 
still keeping g small enough to not lose perturbative control over loops involving 
g.) However the penalty paid in this case is unacceptably large active-sterile mixing 
angles, 9 si. Although wc have not been able to find a choice of parameters within 
this class which satisfies all observational constraints, neither have we been able to 
completely exclude that this is possible. 

3.4 Approximate Symmetries II 

We next turn to a pattern of lepton and flavour symmetries which can produce ac- 
ceptable neutrino phenomenology. We present this model in the spirit of an existence 
proof that such models can be possible within the SLED framework. 

As before we must demand brane-bulk neutrino couplings which conserve a lep- 
ton number in order not to generate too large neutrino masses as we integrate out 
scales between the TeV scale and the sub-eV scale. In this case, however, we assume 
that our couphngs at the TeV scale separately conserve the lepton numbers Lg and 
+ Lt, with couplings 

AW=(0,A,a)" and A!") = . (3.32) 

With this choice all three neutrino masses vanish, and the masslessness of the three 
lightest states is protected by the lepton number invariance. 

We next assume that the renormalization down to the sub-eV scale breaks lepton 
number, but in such a way that the combination Lg — — Lt is less badly broken 
than the other combinations, in which case 

9^0^^ = AW^; = (O, g, gf and = X^-^v = (e, 0, o)" . (3.33) 

As before we take e ~ KniKK/Mg <^ g <^ 1, and so the hghtest nonzero neutrino 
masses are of order egS/r. For simplicity we also assume here the permutation 
symmetry which interchanges the 2nd and 3rd generations, although perturbations 



sin 9 si ~ 



[Oii 



af + 2 
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about this limit are permissible provided they do not move the atmospheric neutrino 
oscillations too far away from maximal mixing. 

With these choices we have the usual massless mass eigenstate 

i/o = ^(o,l,-l,0,---)", (3.34) 

for which Xqj = y^j = sin^so = 0. In this case the unbroken — — Lj. symmetry 
ensures that the two nonzero eigenvalues which vanish as g ^ form a Dirac pair 
with lowest-order mass 

V2egS 



and eigenvectors 



vf ' ' ' ' ' ' ' 



(3.35) 



(3.36) 



to leading order in g. These expressions use the lowest-order result a° = ±a/2 and 



and so 9. 



Xi± = =F^=— and ye± = 

V2q Q 

^o, where as before 



(3.37) 



tan^ 9s ^ g^V , 



(3.38) 



Phenomenologically acceptable masses and mixings are obtained by perturbing 
away from this limit, replacing eqs. (|3.33D with 



9. 



-'2g',g,g 



and 



9a 



ri ' '/I ^ 



/ / 



(3.39) 



where g'/g and e'/e are both 10%, but g'/g — e'/e < 1%. This leads to the lowest-order 
mass matrix 

/ _4e^' eg-2e'g' eg-2e'g'\ 



~{9i-v,-^+9i^9':^ 



S 

r 



(3.40) 



eg - 2e'g' 2ge' 2ge' 
\eg-2e'g' 2ge' 2ge' J 

whose properties we evaluate perturbatively in g' / g and e'/ e. 

With these choices the massless state is not perturbed: 5yUo = and 5vq = 0. On 
the other hand, including the leading and next-to-leading shift in the lightest (but 
nonzero) two perturbed mass eigenvalues gives 



l±V2i--^- 

e 9 



9 



(3.41) 
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where /i^ is given in eq. ( 3.35| ). We work to next-to-leading order in this expression 
since our assumptions imply that the difference g' / g — t' /e is the same order as {g' / gY 
and {e'/ef. 
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Figure 1: The i^e survival probability 
function of L/E in m/MeV, using 
e'/e = 0.1, g'/g = 0.11, g = lO'^, V = 
32, 1/r = 0.01 eV and gP^ = 0.05 eV. 

The corrected values for a± are a± = 



Figure 2: The difference between the 
full survival probability and the survival 
probability computed only using the 'ac- 
tive' states z^o v±. 



^1 



6 



where 



±V2 and 



5 = 2(1 + ^ 

e 9 



(3.42) 



and so the mixing matrix for the lightest three neutrinos becomes to this order 

/c,(-l/V2-5/4) Cs{l/V2-6/A) ^ 

c,{1/2-6/AV2) Cs{1/2 + 6/AV2) 1/^2 . (3.43) 
^c,(l/2- 5/4^2) Cs{l/2 + 6/4V2) -1/V2 j 

Figure 1 shows the survival probability for vacuum oscillations, computed by 
numerically performing the mode sum in the expression 



COS 



Am^L 
2E 



(3.44) 



using a representative set of parameter values and the KK spectrum of a square torus 
(with the zero mode removed). By contrast, Figure 2 shows the contribution of the 
sterile neutrinos to the survival probability. This shows that the oscillations are 
well described by oscillations among the usual three active flavours, with the sterile 
oscillations only introducing an unobservably small 'fine structure'. Furthermore, 
the active mixing angles given by eq. ( 3.43|) can describe the data when g'/g and e'/e 
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are approximately 10% in size, because a 10% shift is sufficient to make the solar 
mixing angle acceptably far from maximal mixing. Yet the ratio of mass splittings 
relevant to solar and atmospheric neutrinos is also acceptable because 

Pf= , ^ -^-KsV2(i-l\^l%. (3.45) 

Because of the factor of 8-\/2 this works provided the difference e'/ e — g' / g is about 
1% of the value of each of e'/e and g' / g separately, which unhke the other choices we 
have made does not seem to follow purely from symmetry grounds. 

Active-sterile mixings are kept acceptable in this model by taking g < 10~^, and 
this in turn requires at the constant K appearing in e = Km.KK/Mg is order 10^ so 
that the mass scales are //± ~ €gS/r ~ Kg/r ~ 1/r ~ 0.01 eV, as required. Since 
K — Mg/Mj, where Mj is the mass scale of the sector which communicates the bulk 
lepton- number breaking to the brane, taking Mg = 10 TeV implies the worrisomely- 
low value Mj ~ 1 GeV. The question whether these choices for K, e' and g' can be 
obtained from a simple extra-dimensional geometry is presently under study. 



4. Conclusions 

The coincidence between the values of the recently-discovered neutrino masses and 
the Dark Energy density suggests the possibility of a deep connection between these 
quantities — a possibility which can be more sharply posed within the SLED pro- 
posal for resolving the cosmological constant problem. Embedding extra- dimensional 
neutrino models into SLED opens up many attractive features, largely because the 
supersymmetry of the bulk physics necessarily constrains the otherwise arbitrary 
properties of the bulk sterile neutrinos. In particular, it removes the need for the 
ad hoc assumption of vanishing bulk fermion masses, since supersymmetry requires 
these masses to vanish. As is also true for non-supersymmetric models, brane-bulk 
neutrino mixing provides a natural extension of lepton number to bulk fermions, by 
embedding it within local Lorentz rotations of the large extra dimensions. 

In this paper we provide an example of a model which realizes this connection in 
a concrete way, thereby motivating a more detailed examination of neutrino models 
within the SLED context. The possibihty of so doing comes as a surprise, since 
past experience of neutrino models involving large extra dimensions were driven 
by sterile-neutrino constraints towards extra dimensions having properties which are 
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inconsistent with the requirements of the SLED picture for the cosmological constant 
problem. In particular, constraints on catastrophic neutrino emission into the bulk 
typically drive model builders to choose only one dimension to exist near the eV 
scale, since otherwise the higher-dimensional phase space can make this emission 
rate much too large. 

By contrast, we here find that astrophysical and cosmological bulk-neutrino emis- 
sion can be acceptably small provided that the brane-bulk mixing is itself small 
(51/277 < 10~^). The somewhat surprising result that this is possible follows largely 
from the observation that back-reaction of the stress-energy of the brane typically 
removes the massless fermion KK mode. (For example, many co-dimension 2 branes 
induce a conical singularity in the extra-dimensional geometry at the brane position 
and the boundary conditions which this conical defect imposes typically removes the 
massless bulk fermion modes.) The removal of these zero modes is crucial for suc- 
cessful neutrino phenomenology, since even a small brane-bulk couplings typically 
induces large mixings amongst any degenerate massless bulk and brane states. If 
no bulk zero modes exist then active-sterile neutrino mixings remain small for small 
brane-bulk couplings, and the large mixings amongst the active neutrinos can explain 
the observed neutrino oscillation pattern. 

We construct a model having acceptable phenomenology by assuming a relatively 
mild flavour-dependence for the brane-bulk couplings, and by taking these couplings 
to be sufficiently small. Technically natural small neutrino masses are obtained, with 
approximate lepton number protecting the small masses of the lightest neutrinos. 
This is possible because lepton number extends into the extra dimensions as a local 
Lorentz transformation, and for generic internal geometries it is naturally broken 
at the Kaluza-Klein scale. A potentially unpleasant feature of the small brane-bulk 
couplings which we assume is that the amount of lepton number breaking which is 
necessary to get acceptable neutrino oscillations appears to require a relatively large 
coupling of bulk lepton-number breaking to the brane neutrinos. 

Although it is not yet clear whether such couplings can be generated from real 
compactifications, we believe the success of the models described herein to motivate 
more detailed studies of the possibilities offered by neutrino model building within 
the SLED framework. 
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A. Diagonalization for Large Mixings 

In this appendix we present the exact results for the eigenvectors and eigenvalues for 
the mass matrix in the two cases discussed in the main text. 

A.l Approximate Symmetries I 



We take here the symmetry choice of Section O and impose = for simplicity, 
implying 

^(+) = (ei, 62, 62)^ ^("^ = {g, g, gf 

The eigenvectors and eigenstates for the mass matrix in this case come in two distinct 
kinds: 

• There is an eigenstate whose eigenvalue is exactly zero: 

1/0 = ;^ (0,1, -1,0, 0,0,0,...) 

• All other eigenstates have nonzero eigenvalues given by fii = \zi\/r, where the 
quantity Zi is the solution of the following transcendental equation 



-2gS',e^ +[1~V: {g' + ej)] z?) [-4gS:e2 + [l - 2V[ + e^)] 

-2 [gS[{t^ + £2) + V[{g^ + £162)^.) ' = , (A.l) 

where 

S[ = Y,Su with S[,^^^ and V[ = Y,Pk with P'h^^^ 

(A.2) 



_ ^ 2 



The corresponding eigenvectors are 

z/j = Ni (aj, 1, 1, xh, uii, X2i, y2i, xsi, y^u ■■■) , (A.3) 
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where 

2gS[{ti + t2) + 2V[z,{g^ + eie2) 



a. 



(1 - {g' + eD^D - 2ge,S[ ' 
xti = gS[i (ttj + 2) + ZiP'^i {2t2 + ttiei) , 

yu = (262 + caei) + ^z^P;, {a^ + 2) (A.4) 
and the overall normalization factor is given in terms of these by 

= I ^ . (A.5) 

The solutions to the eigenvalue equation ( |A.1|) come in the following two different 
types: 

• Two light eigenvalues that in the limit of z± <C c\ are: 

1/2- 



z± = -gS- 

where 



(ei + 262) ± ((ei + 262)2 + 2(6i - 62)'n, 



(A.6) 



= 1 + {?,g^ + el + 2el)V + 2(6i - e^fg^V^ . 
The corresponding eigenvectors are obtained by using eq. (|A.6|) in eq. (|A.3|) . 

• An infinite tower of eigenvalues, = \zi±\/r which come in pairs such that 
Zii± —>■ ±ci as g —>■ 0. This corresponds to a full tower of Kaluza Klein states 
whose masses are slightly perturbed away from the zeroeth-order result \ce\/r 
by the brane-bulk couplings. 

A. 2 Approximate Symmetries II 

We now present expressions for the eigenvalues and eigenvectors of the mass matrix 



in the case discussed in Section 3.4. It is useful here to distinguish the following two 



subcases, as discussed in the main text: 

a) Unperturbed case: g'^ = {0,g,g)'^ and g~ = (6,0,0)"^ 



b) Perturbed case: g'^ = {—2g\g,gY and g = {e,e',e'Y 
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Case a: Unbroken Lepton Symmetry 

The eigenvectors in this case resemble those found above: 

• There is an eigenstate with zero eigenvalue: 

1/0 = ^(0,1,-1,0,0,0,0,...) . 

• The rest of the eigenvectors have the form: 

Vi = Ni (ttj, 1, 1, xii, yii, X2i, y2i, x^i, ysi, ...) 

where 

2Sige 

OLi 



xii = S'^^aie + ZiPl,i2g, 

yii = S[i2g + ZiP'^^aie (A. 7) 

and the normalization factor A^j is defined as in eq.( |A.5| ). As before, Zi is related 
to the mass eigenvalue /Xj by /ij = \zi\/r, with Zi defined as the solutions to the 
eigenvalue equation: 

(1 - (1 - - 2 S['gh' = , (A.8) 



with V' and S' defined as in eq. ( [A.2|) . The solutions to this equation come in two 
types: 

• Two light eigenvalues given by 

V 

with 

Ui = -2- 2{e^ + 2g^)V - Ae'^g^V^ 

p. = ~{egfS' (A.9) 

where we have simplified the result as appropriate for those eigenvalues satisfy- 
ing zf <C c|. The corresponding neutrino masses are in this case do not depend 
on the label ± since they have the common mass /x± = |-z±|/t. 

• Also in this case we have the tower of eigenstates whose eigenvalues (for small 
g) are perturbations to zi± ~ ±q. 
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Case b: Broken Lepton Symmetry 

Perturbing the above results with the addition of e' and g' leads to a phenomeno- 
logically viable model. In this case we get the same structure of eigenvectors and 
eigenvalues, with the important difference that now the degeneracy between the 
lightest massive neutrinos is broken: 

• One eigenstate remains massless, with eigenvector 

z/o = ^(0,1,-1,0,0,0,0,...). 

• The rest of eigenvectors are given by: 

Vi = Ni {ai, 1, 1, xii, yu, X2i, y2i, xsi, y^i, ...) , (A.IO) 

where 

2S[ {ge - 2g'e') + 2V[z, {-2gg' + ee') 

fy ■ = 

z,{l-{Ag'^ + e^)V'^+Ag'eS[ ' 

xti = 5; (a^e + 2e') + z,P;, {2g - 2aig') , 

ya = {2g - 2aig') + z,P',, {ait + 2e') (A. 11) 

with the normalization factor defined as in eq. ( |A.5[) . The eigenvalue equation in 
this case reads: 

[z, [I - {Ag'^ + e^)^:] + Ag'eS') {z, [l - 2(/ + e'^)p;] - Age's) 

-2{S[ {gt - 2g't') + V[z, {-2gg' + ee'))' = . (A. 12) 

The solutions to this equation have the same structure as before. 

• There are two light non-degenerated eigenvalues given by 

_Vi± a/ vf + AuiPi 



(A.13) 



with 



Vi = - A{eg' - e'g)S 

Ui = -2- 2(e2 + 2(e'2 + g^ + 2g'^))V - A{eg + 2eV)'^' 

= -{eg + 2t'g'fS\ (A.14) 

where we simplify by assuming zf <^ c^. As usual, the corresponding neutrino 
mass is /i± = l^^il/r. The eigenvectors are obtained by using eq. (|A.13 | ) in 
eq.(CT). 
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• There is a tower of eigenstates organized in pairs, having eigenvalues which are 
perturbations of Zi± ~ ±q. 
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